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Abstract 


In  sarin  (isopropyl  methylphosphonofluoridate)  there  are  ten  nonequivalent  CH  oscil¬ 
lators.  Ab  initio  calculations  at  the  HF  /  6-31I++G(2d32p)  level  have  been  used  to 
determine  the  vapour  phase  local  mode  parameters,  u;  and  lux,  for  each  oscillator  in 
the  two  spectrally  significant  conformers  of  sarin,  as  well  as  inter-oscillator  coupling 
parameters.  These  above  parameters,  in  conjunction  with  dipole  moment  functions 
derived  from  ab  initio  calculations,  were  used  to  perform  harmonically  coupled  anhar- 
monic  oscillator  (HCAO)  calculations,  thereby  enabling  the  simulation  of  vibrational 
overtone  spectral  regions  in  a  room- temperature  sample  of  sarin.  It  was  determined 
that  the  eomputationally-intensive  HCAO  approach  is  necessary  to  predict  the  lower 
vibrational  overtone  regions  (i,  e.,  first  to  third  overtones)  as  a  simpler  “non-HCAO” 
approach  (which  does  not  allow  pairwise  harmonic  coupling  among  adjacent  oscilla¬ 
tors)  failed  to  accurately  reproduce  the  HCAO-simulated  spectral  regions. 

The  present  work,  which  was  carried  out  without  recourse  to  the  experimental  sarin 
spectral  regions,  illustrates  that  it  is  currently  feasible  to  predict  the  absorption 
spectra  of  species  which  are  difficult  to  synthesize,  handle,  or  otherwise  acquire.  In 
addition  to  their  utility  in  guiding  experimental  investigations,  the  simulated  overtone 
spectral  regions  will  be  necessary  to  correctly  assign  experimental  overtone  spectra 
owing  to  the  large  number  of  similar  but  nonequivalent  CH  oscillators  present  in 
sarin. 

Resume 


II  y  a  dix  oscillateurs  CH  non  equivalents  dans  le  sarin  (metbylfluorophosphonate 
d’isopropyle).  Les  calculs  Ab  initio  au  niveau  HF  /  6-311-t-+G(2dJ2p)  ont  ete  utilises 
pour  determiner  les  parametres  en  mode  local  de  la  phase  gazeuse,  et  wx,  pour 
chaque  oscillateur  dans  les  deux  conformeres  de  sarin  specifiquement  signifiants  spec- 
tralement,  ainsi  que  les  parametres  de  couplage  inter-oscillateur s .  Les  parametres 
ci-dessus,  en  conjunction  avec  les  fonctions  moment  dipole,  derivees  des  calculs  ab 
initio ,  ont  ete  utilises  pour  effectuer  des  calculs  harmoniquement  couples  d 'oscilla¬ 
teurs  non  harmoniques  (HCAO),  permettant  ainsi  la  simulation  de  regions  spectrales 
ayant  une  resonance  vibrationnelle  dans  un  echantillon  de  sarin  de  temperature  am- 
biante.  On  a  determine  que  la  methode  HCAO  aux  calculs  intensifs  est  necessaire 
pour  predire  les  regions  ayant  une  resonnance  vibrationnelle  plus  basse  (par  ex.  :  de 
la  premiere  a  la  troisieme  resonnance)  parce  que  la  methode  <Cnon  HCAO»  plus 
simple  (qui  ne  permet  pas  le  couplage  hannonique  en  paires  par  mi  les  oscillateurs 
adjacents)  n?a  pas  reussi  a  reproduce  avec  exactitude  les  regions  spectrales  simulees 
par  la  HCAO. 
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Les  travaux  presents,  qui  ont  ete  effect  ues  sans  avoir  recours  aux  regions  spec t rales 
de  sarin  experimental,  illus t rent  qu’il  est  actuellement  faisable  de  predire  le  spectre 
d 'absorption  des  especes  qui  sont  difficiles  a  svnthetiser,  manipuler  on  encore  acquerir. 
En  plus  de  leur  utilite  a  guider  les  investigations  experiment  ales,  les  regions  spec- 
t rales  de  resonnance  simulees  seront  necessaries  pour  assigner  correctement  le  spectre 
experimental  des  resonnanees  d  a  la  grande  quantite  d'oscillateurs  CH  non  equivalents 
presents  dans  le  sarin. 
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Executive  summary 


Computational  Simulation  of  Vibrational  Overtone 
Spectral  Regions:  Sarin 

Michael  W.  P.  Petryk;  DRDC  Suffield  TR  2006-220;  Defence  R&D  Canada  - 

Suffield;  December  2006. 

Background:  There  is  a  strong  interest  in  developing  the  ability  to  detect  chemical 
threats  at  a  standoff,  that  is,  beyond  the  effective  range  of  the  hazards  which  these 
agents  pose.  The  ability  to  safely  interrogate  a  suspect  contaminant  at  a  distance 
implies  the  use  of  some  form  of  spectroscopy,  in  which  light  interacts  with  the  threat 
agent  and  brings  about  a  change  (such  as  a  spectral  absorption  feature)  which  can 
be  detected  at  a  distance.  Historically,  efforts  in  the  standoff  detection  of  chemical 
warfare  agents  (CWAs)  have  focussed  on  the  use  of  infrared  (IR)  light  in  the  2.5-25 
wavelength  region.  However,  recent  advances  the  telecommunications  industry  such 
as  the  new  availability  of  small,  inexpensive,  and  efficient  near-infrared  (NIR)  diode 
lasers  now  encourage  further  investigations  into  the  possibility  of  detecting  the  CWAs 
spectroscopically  in  the  NIR  (i.  e.,  0.8-25  /.im)  and  visible  regions  (400-750  nm). 

Principal  results:  It  is  often  desirable  to  be  able  to  predict,  without  recourse  to  any 
experimental  data,  the  spectral  absorption  features  of  CWAs.  This  is  especially  im¬ 
portant  if  the  CWA  in  question  is  difficult  to  handle,  synthesize,  or  otherwise  obtain. 
In  this  report  it  is  demonstrated  that  it  is  feasible  to  predict  the  spectral  absorp¬ 
tion  features  of  sarin  (isopropyl  methylphosphonofluoridate)  by  relying  only  upon 
computer  models,  without  any  recourse  to  experimental  data.  Two  computational 
models  have  been  evaluated  in  this  report,  only  one  of  which  was  found  to  produce 
satisfactory  results. 

Significance  of  results:  This  report  illustrates  that  it  is  feasible  to  predict  (without 
recourse  to  the  experimental  data)  the  absorption  spectra  of  species  which  are  difficult 
to  synthesize,  handle,  or  otherwise  obtain.  In  addition  to  their  utility  in  guiding 
experimental  investigations  and  equipment  design,  the  simulated  spectral  regions 
presented  herein  will  be  necessary  to  correctly  assign  and  interpret  transitions  in  the 
experimental  spectral  regions. 

Future  work:  Further  validation  of  the  approach  to  spectral  simulation  reported 
herein  should  be  undertaken  by  acquiring  the  experimental  spectral  regions  of  several 
CWAs  in  the  vapour  and  liquid  phases. 
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The  ab  initio  calculations  which  were  used  in  the  research  presented  herein  were  car¬ 
ried  out  on  a  publicly- funded  academic  computer  network.  It  would  be  advantageous 
for  Defence  Research  and  Development  -  Suffield  to  acquire  the  necessary  hardware 
and  software  in  order  to  be  able  to  carry  out  such  research  in  a  secure  environment. 
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Contexte  :  II  y  existe  un  vif  interet  a  mettre  au  point  la  capacite  a  detecter  des  me¬ 
naces  chimiques  a  distance,  c’est  a  dire  au-dela  de  la  portee  efficace  des  menaces  que 
posent  ces  agents.  La  capacite  a  interroger  a  distance  un  contaminant  suspect  en  toute 
securite  implique  l’utilisation  dime  forme  de  spectroscopie  dans  laquelle  la  lumiere  in- 
teragit  avec  1’ agent  de  menace  et  apporte  un  changement  (tel  qu'une  caracteristique 
d’absorption  spectrale)  pouvant  etre  detectee  a  distance.  Les  efforts  concernant  la 
detection  a  distance  d’agents  de  guerre  chimiques  ont  historiquement  focalise  sur  l’uti- 
lisation  de  la  lumiere  infrarouge  (IR)  d’une  longueur  d’onde  d'environ  2,5  a  25  pm. 
Les  progres  recents  dans  le  domaine  de  l’industrie  des  telecommunications  tels  que  la 
nouvelle  disponibilite  de  petits  lasers  a  diode  proche  infrarouge  (NIR)  economiques 
et  efficaces  encouragent  cependant  maintenant  a  examiner  plus  profondement  la  pos¬ 
sibility  de  detecter  les  armes  de  guerre  chimiques  spectroscopiquement  dans  le  NIR 
(par  ex.  :  de  0,8  a  25  pm)  et  dans  des  regions  visibles  (de  400  a  750  nm). 

Resultats  principaux  :  II  est  souvent  desirable  d’etre  en  mesure  de  predire  les 
caracteristiques  d’absorption  spectrale  d’une  arme  chimique  sans  avoir  recours  a  des 
donnees  experimentales.  Ceci  est  surtout  important  si  l’arme  chimique  en  question  est 
difficile  a  manipuler,  synthetiser  ou  obtenir.  Ce  rapport  demontre  qu’il  est  faisable  de 
predire  les  caracteristique  d’absorption  spectrale  du  sarin  (methylfluorophosphonate 
d’isopropyle)  en  se  fiant  seulement  a  des  modeles  informatises  et  sans  aucun  recours  a 
des  donnees  experimentales.  Deux  modeles  de  calcul  ont  ete  evalues  dans  ce  rapport 
mais  un  seul  a  produit  des  resultats  satisfaisants. 

Portee  des  resultats  :  Ce  rapport  illustre  qu’il  est  faisable  de  predire  le  spectre 
d’absorption  des  especes  difhciles  a  synthetiser,  manipuler  ou  obtenir  (sans  avoir  re¬ 
cours  a  des  donnees  experimentales).  En  plus  de  leur  utilite  a  guider  les  investigations 
experimentales  et  la  conception  de  l’equipement,  les  regions  spectrales  simulees  men- 
tionnees  seront  necessaires  pour  assigner  et  interpreter  correctement  les  transitions 
dans  les  regions  spectrales  experimentales. 

Travaux  futurs  :  La  validation  future  de  la  methode  de  simulation  spectrale  docu- 
mentee  ici  devrait  etre  entreprise  en  acquerant  les  regions  spectrales  experimentales 
de  plusieurs  armes  chimiques  en  phase  gazeuse  et  liquide. 
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Les  calculs  ah  initio  utilises  dans  la  recherche  presentee  ici  ont  ete  effect  ues  sur  un 
reseau  informatique  pedagogique  finance  par  le  gouvernement.  II  serait  avantageux, 
pour  Recherche  et  developpement  pour  la  defense  -  Suffield,  d'acquerir  le  materiel  et 
les  logiciels  requis  pour  effectuer  une  telle  recherche  dans  un  milieu  securitaire. 
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1  Introduction 


Vibrational  stretching  overtone  transitions  of  CH  bonds  (which  occur  approximately 
at  frequencies  of  n  x  3000  cm-1  where  n  is  an  integer  which  is  greater  than  or  equal  to 
two)  are  potentially  exploitable  for  the  portable  spectroscopic  detection  of  chemical 
warfare  agents  (CWAs)  as  many  such  transitions  occur  in  the  near-infrared  (NIR)  and 
visible  regions  of  the  spectrum  where  small,  efficient  diode  lasers  operate.  Further, 
uncooled  detectors  which  operate  in  the  NIR  and  visible  regions  are  far  more  sensitive 
than  those  which  operate  in  the  more  commonly  exploited  fundamental  region  (t.  e., 
ca.  450-3500  cm-1).  For  example,  uncooled  indium  gallium  arsenide  (InGaAs)  diode 
detectors  operate  in  the  5800-12500  cm-1  region  and  are  approximately  six  hundred 
times  more  sensitive  than  cryogenically  cooled,  wide-band  mercury  cadmium  telluride 
(MOT)  detectors  which  are  commonly  used  in  fundamental  vibrational  spectroscopy 
and  operate  in  the  420-10000  cm-1  region.  If  cryogenic  cooling  is  impractical  or  im¬ 
possible,  uncooled  deuterated  triglycine  sulfate  (DTGS)  detectors  are  available  which 
operate  in  approximately  the  same  region  as  cooled  MCT  detectors.  However,  the 
sensitivities  of  uncooled  DTGS  detectors  are  approximately  twenty-five  times  lower 
than  those  of  MCT  detectors.  Unfortunately,  vibrational  overtone  transition  inten¬ 
sities  decrease  rapidly  with  increasing  vibrational  excitation  such  that  enhancements 
realized  by  increases  in  portable  laser  power  and  detector  sensitivities  in  the  lower 
overtone  regions  are  negated  at  CH  overtones  which  lie  approximately  at  or  above  the 
fourth  overtone.  This  report  deals  with  the  computational  modelling  of  vibrational 
overtone  transitions  in  the  nerve  agent  sarin  (isopropyl  methylphosphonofiuoridate). 
No  attempt  is  made  to  propose  specific  spectroscopic  techniques  for  the  detection  of 
vibrational  overtone  transitions  in  CWAs  as  such  specific  techniques  will  have  to  be 
chosen  based  upon  the  intended  application  (e,  g.,  point  detection,  standoff  detection 
suitable  surveying  large  areas,  etc.). 

The  vibrational  overtone  stretching  transitions  of  XH  bonds  (where  X  =  C,  N,  O, 
etc.)  are  dominated  by  transitions  into  states  where  all  of  their  vibrational  energy 
is  localized  within  one  of  a  set  of  equivalent  XH  bonds.  The  energy  localized  within 
these  local  modes  of  vibration  [1-5]  makes  overtone  spectroscopy  a  very  sensitive 
probe  for  detecting  small  changes  in  the  properties  of  XH  bonds,  such  as  molecular 
conformations  [6,7]  and  bond  lengths  [8-10], 

The  molecular  Hamiltonian  can  be  expressed  in  a  local  mode  basis.  Such  a  repre¬ 
sentation  leads  to  XH  stretch  terms  where  off-diagonal  coupling  elements  are  small 
compared  to  the  (anharmonic)  diagonal  cubic  and  quartic  terms  [11].  The  contribu¬ 
tions  of  the  diagonal  terms  to  the  molecular  Hamiltonian  become  increasingly  dom¬ 
inant  with  increasing  vibrational  excitation.  As  a  result,  it  is  possible  to  treat  the 
stretching  potential  of  a  molecule  as  a  collection  of  isolated  anharmonic  oscillators 
in  which  inter -oscillator  coupling  is  treated  as  a  perturbation  [12].  Often  such  an- 
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harmonic  oscillators  are  approximated  using  Morse  oscillators  since  the  energies  and 
matrix  elements  of  the  wavef  unctions  associated  with  the  latter  have  been  solved 
analytically  [4,13]. 

If  the  stretching  potential  of  a  molecule  is  approximated  as  a  collection  of  isolated 
(Morse)  anharmonic  oscillators,  the  basis  states  of  a  molecule  comprising  n  A'H  Morse 
oscillators  can  be  written  as  the  product  of  individual  one  dimensional  Morse  wave- 
functions  \vi)\v2)  -  *  *  \Vi)  -  ■  ■  Wn)  or  simply  as  \vuv2 , . . . ,  v*, . . . ,  vn },  where  v{  is  the 
number  of  vibrational  quanta  in  the  zlh  oscillator.  Those  states  which  carry  intensity 
from  the  v  —  0  ground  vibrational  state  (i.  e.,  spectrally  “bright1'  states),  are  almost 
exclusively  pure  local  mode  states  which  are  described  by  a  linear  combination  of 
components  in  which  all  vibrational  energy  is  localized  within  one  oscillator  (e.  g., 
|tf,  0, . . .  ,0),  10,1/, ...  *0),  etc.).  In  the  case  where  the  local  methyl  group  has  C\ 
symmetry  (e.  sarin),  the  pure  local  mode  wavefunctions  of  the  three  uncoupled 
methyl  XH  oscillators  are  denoted  |^i)|0)|0),  |Q))Vj)|0),  and  |0))0}|i'jt).  In  the  case 
of  a  CH  oscillator,  a  transition  out  of  the  ground  state  v  —  0  into  a  state  excited 
by  n  vibrational  quanta  is  often  referred  to  by  the  change  in  quanta  and  denoted 
Avqu  ~  n -  The  local  mode  description  has  been  extended  in  the  harmonically  cou¬ 
pled  anharmonic  oscillator  (HCAO)  model  [5, 14—17]  wherein  molecules  are  treated 
as  a  collection  of  Morse  oscillators  which  are  harmonically  coupled. 

The  HCAO  model  can  accurately  predict  the  intensities  of  light  absorption  which 
accompany  vibrational  overtone  transitions,  both  in  terms  of  the  relative  absorbance 
intensities  of  oscillators  within  a  molecule  [18-21],  and  in  terms  of  the  absolute  ab¬ 
sorbance  intensities  [22-28].  The  HCAO  model,  used  in  conjunction  with  ab  initio 
energy  and  dipole  moment  calculations,  has  allowed  the  accurate  prediction  of  over¬ 
tone  spectra  for  species  which  have  not  been  observed  experimentally  [29,30].  The 
transition  from  normal  mode  vibrational  behaviour  to  local  mode  vibrational  be¬ 
haviour  (the  HCAO  model  is  based  upon  the  latter)  in  XH  stretching  oscillators 
occurs  at  approximately  Auch  =  3.  Simulated  CH  vibrational  overtone  spectra] 
regions  will  be  presented  for  sarin  in  the  AuCh  =  3  —  6  regime. 

In  sarin  there  are  ten  nonequivalent  CH  oscillators.  In  this  report  ab  initio  calcula¬ 
tions  will  be  used  to  determine  the  local  mode  parameters,  u}  and  ujx,  for  each  oscilla¬ 
tor  in  all  spectrally  significant  conformers  of  sarin,  as  well  as  inter-oscillator  coupling 
parameters.  These  above  parameters,  in  conjunction  with  dipole  moment  functions 
derived  from  ab  initio  calculations,  will  be  used  to  perform  HCAO  calculations,  en¬ 
abling  the  simulation  of  vibrational  overtone  spectral  regions  in  a  room- temperature 
sample  of  sarin. 
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2  Theory  and  calculations 

2.1  Anharmonic  oscillators  and  local  modes  of 
vibration 


Vibrational  overtone  spectra  arising  from  CH  stretching  transitions  can  be  interpreted 
within  the  local  mode  model  of  molecular  vibration  (2,3].  The  CH  stretching  potential 
is  often  approximated  by  a  Morse  potential  [12] 

V{q)  =  Dt{  l-e““’)2  (1) 


where  De  is  the  depth  of  the  potential  energy  well,  o  is  the  Morse  scaling  factor,  and 
q  is  displacement  from  the  equilibrium  CH  bond  length.  The  depth  of  the  potential 
energy  well  is  given  by 

Dt-  (2) 


4 ujx 


where  the  terms  u  and  u5x  are,  respectively,  the  pure  local  mode  frequency  and 
anharmonicity  (see  below).  The  Morse  scaling  factor  o  is  given  by 


(3) 


where  /i  is  the  reduced  mass  of  the  CH  oscillator.  The  energy  eigenvalues,  £V1  of  a 
Morse  oscillator  can  be  expressed  (in  cm-1)  by  the  expression  [12,31] 


By  =  u>  Tv  +  ^  -  Dx  +  0 


(4) 


This  equation  can  be  rearranged  to  give  the  frequency  of  an  overtone  transition  out 
of  the  ground  state,  Swo*  in  the  form  [12] 

uv^o  =  vui  —  (v2  +  v)u)x  .  (5) 


Experimental  values  of  the  Morse  parameters  w  and  £x  can  be  obtained  either  from 
a  Birge-Sponer  plot  of  vs.  v  or  from  a  quadratic  fit  of  vs.  v.  Ab  initio 

calculations  can  also  be  used  to  predict  local  mode  parameters  [32,33]. 


An  anharmonic  oscillator  model  for  molecular  vibrations  which  neglects  coupling 
between  adjacent  oscillators  (as  opposed  to  the  HCAO  approach,  see  Section  2.2), 
hereinafter  referred  to  as  the  “non- HCAO”  approach,  can  be  used  to  predict  the 
frequencies  of  vibrational  overtone  transitions  using  ab  initio-calculated  potential 
energy  surfaces  (PESs)  and  equations  1-3  and  5.  Such  an  approach  is  known  to 
produce  systematic  errors,  overestimating  transition  frequencies  by  approximately 
10-30  cm-1  for  lower  overtones  owing  to  the  neglect  of  interoscillator  coupling  (see 
Section  3.3  for  more  details).  By  way  of  comparison,  the  magnitude  of  the  error  in  a 
typical  HC AO-calculated  transition  frequency  is  approximately  1-3  cm'1. 
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2.2  Coupled  methyl  oscillators 

The  frequencies  of  transitions  out  of  the  ground  |0)  state  to  the  excited  state  \v) 
for  an  isolated  CH  oscillator  are  given  by  equation  5,  In  the  case  of  coupled  CH 
oscillators  on  a  methyl  group  with  local  Cj  symmetry  (e.  sarin),  transitions  out  of 
the  ground  |0)|0)|0)  state  into  the  states  |^}|0)|0),  |0)|y/)|0),  or  |0)|t)})r;jt)  have  an 
equivalent  energy  expression,  given  in  terms  of  the  molecular  Hamiltonian  H,  which 
is  [5] 

-  -  f;,l°>l°>  =  +«,)o*, 

t  t 

”  Y  Ty  (a,7 aj  +  *T a~ )  (6) 

where  E|0)|o)|n)  is  the  energy  of  the  vibrational  ground  state,  a+  and  a~  are  (respec¬ 
tively)  tlie  creation  and  annihilation  operators,  and  7J  .  is  the  intramanifold  coupling 
parameter  [5] 

iii  =  {la  -  <t>a)  y/SiUj  (7) 

where  jl3  and  4>ij  are,  respectively,  the  kinetic  and  potential  energy  coupling  terms. 
These  terms  can  be  expressed  as  [5] 


,  1  Fij 

V'j  -  2F«Fjj 

(8) 

1  Gy 

7,J  2  GaGjj 

(9) 

where  F  and  G  represent  the  Wilson  matrix  elements  [34].  The  diagonal  force  con¬ 
stants  Fu  and  F}]  can  be  derived  from  ab  initio  PESs  along  single  stretch  coordinates 
while  the  non-diagonal  force  constant  F^  can  be  obtained  from  ab  initio  PESs  along 
two  stretch  coordinates  [35],  The  kinetic  energy  coupling  term  7^  is  a  function  of 
molecular  conformation  only  and  can  be  expressed  as  [5] 

7,i  =-^c°s(0)  0  +  ~)  <10> 

where  me  and  mn  denote,  respectively,  the  atomic  masses  of  C  and  H,  The  value  of 
the  H-C-H  bond  angle  0  if  often  calculated  ab  initio. 

2.3  Oscillator  strength 

Experimental  vibrational  overtone  transition  intensities  can  he  expressed  in  terms  of 
the  dimensionless  oscillator  strength.  /osc  [36].  By  applying  the  ideal  gas  law  and 
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using  appropriate  physical  constants  [37]  an  expression  for  can  be  obtained  in 
the  form 

fosc  =  2,6935  x  10-9(Torr  m  cm)  —  J  A{vtg)A{yeg)  (11) 

where  T  is  the  gas  temperature  {in  Kelvin),  p  is  the  pressure  of  the  gas  (in  Torr),  l  is 
the  path  length  (in  meters),  and  A(veg)  is  the  absorbance  as  a  function  of  frequency, 
i>tg,  expressed  in  wavenumbers  (cm'1)  for  a  transition  from  ground  state  g  to  excited 
state  e.  The  integration  is  carried  out  over  the  frequency  range  of  the  absorption 
feature. 

The  oscillator  strength  of  an  overtone  transition  can  be  calculated  from  the  expres¬ 
sion  [36] 

/ce  =  4.70175  x  10'7(cm  D-a)t>efl|/ies|2  (12) 

where  the  term  4.70175  x  10-7  denotes  a  collection  of  physical  constants  (taken  from 
Reference  [37]),  ueg  denotes  the  transition  frequency  (in  cm-1),  and  peg  is  the  tran¬ 
sition  dipole  moment  matrix  element.  The  transition  dipole  moment  matrix  element 
is  given  by 

=  <e|£]  g)  (13) 

where  e  and  g  are,  respectively,  the  excited  and  ground  state  wavefunctions  and  fl 
denotes  the  dipole  moment  function  (DMF)  expressed  in  Debye  (see  equation  14, 
below).  The  wavefunction  of  a  Morse  oscillator  excited  by  v  vibrational  quanta  can 
be  calculated  from  the  local  mode  parameters  (5  and  ux )  and  the  reduced  mass  of 
the  oscillator,  p  [12]. 

2.4  The  dipole  moment  function 

Within  the  HCAO  model,  coupling  between  adjacent  CH  oscillators  is  treated  as 
a  pairwise  interaction  between  the  oscillator  CHa  along  internal  displacement  coor¬ 
dinate  qB  and  oscillator  CH*,  along  internal  displacement  coordinate  g&.  Ah  initio 
calculations  are  used  [18.19,38-40]  to  obtain  discrete  dipole  moment  values  as  func¬ 
tions  of  qa  and  q t,  In  this  report  the  dipole  moment  values  are  calculated  as  qa  and  q ^ 
are  compressed  or  stretched  about  their  equilibrium  CH  bond  lengths  rc  by  ±  0.30  A 
in  increments  of  0,05  A  with  all  other  molecular  coordinates  constrained  to  their 
equilibrium  values.1  The  continuous  DMF  of  two  coupled  oscillators  along  internal 
stretch  coordinates  qa  and  qt  is  approximated  by  a  Taylor  series  expansion  taken  up 
to  sixth  order 

£(9«>9fc)  =  AirfWb  (14) 

<J=1 

'Examples  of  z- mat  rices  which  correspond  to  the  three  lowest  conformers  of  sarin  are  presented 
in  Annex  A. 
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where  /2y  are  the  DMF  derivatives 


1 

i'-j'-  dqidqi 


(15) 


obtained  from  a  fit  of  ab  initio-calculated  dipole  moment  values  to  a  polynomial  (the 
subscript  eg  denotes  that  all  other  molecular  coordinates  are  constrained  to  their 
equilibrium  values).  Terms  in  equation  14  include  the  pairwise  mixed  derivatives 
d2fl/dqadqt,t  d^fl/dqadq2,  and  d3jl/dq2dqt>.  It  is  important  to  include  terms  beyond 
those  which  are  linear  in  g  since,  for  example,  the  second  order  derivatives  d2jl/dql 
and  d2fi/dq 2  contribute  significantly  more  to  overtone  transition  intensities  than  do 
the  linear  terms.  It  is  known  that  “non-HCAO”  intensities  (where  inter-oscillator 
coupling  is  neglected)  are  systematically  incorrect  (5,14).  The  magnitude  of  the  “non- 
HCAO”  intensity  errors  are  such  that  “non-HCAO”  intensities  are  not  presented  in 
this  report. 


It  is  known  that  HCAO  calculations  better  reproduce  experimental  overtone  transi¬ 
tion  intensities  if  they  use  ab  initio  DMFs  where  computer  resources  have  been  spent 
on  larger  basis  sets  rather  than  on  increased  levels  of  electron  correlation  [41],  espe¬ 
cially  if  the  molecule  does  not  contain  any  unsaturated  bonds.  For  transitions  to  the 
second  vibrational  overtone  and  higher  electron  correlation  effects  on  the  DMF  are 
negligible  [22,41,42].  In  general,  it  has  been  found  that  Hartree-Fock  (HF)  theory 
gives  excellent  DMFs  for  the  purpose  of  HCAO  calculations  of  CH  stretching  overtone 
intensities  [22,41], 


2.5  Computational  details 

.46  initio  calculations  were  carried  out  using  Gaussian  03,  Revision  CX)2  [43]  on  an  HP 
AlphaServer  ES40  cluster  with  830  MHz  and  670  MHz  64-bit  processors.  Calculations 
used  all  Gaussian  defaults  save  that  Gaussian  overlay  option  IOP( 3/32=2)  was  used 
with  all  calculations  to  prevent  the  possible  reduction  of  the  expansion  set.  All 
calculations  were  carried  out  using  SCF=Tight  as  the  6-31l4H-G(2d,2p)  basis  set 
used  in  this  Teport  contains  diffuse  functions. 

For  each  confer mer  of  interest  single  point  energy  calculations  were  carried  out  as 
functions  of  the  GH  stretch  coordinate  g  in  the  displacement  range  q  €  [—0.30, 0.30]  A 
in  steps  of  0.05  A  (a  total  of  13  calculations  per  unpaired  stretch  coordinate).  In 
sarin  there  is  one  “isolated”  (unpaired)  CH  stretch  coordinate  and  three  {coupled) 
CH  stretch  coordinates  on  each  of  three  methyl  groups.  Within  the  HCAO  model, 
the  methyl  CH  oscillators  give  rise  to  a  total  of  nine  (asymmetric)  pairwise  CH 
stretch  interactions.  Thus,  for  each  conformer  a  total  ofI3xl3x9+13  —  1534 
ab  initio  calculations  were  carried  out.  The  total  computation  time  required  for 
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1534  single  point  (energy  and  dipole  moment)  energy  calculations  on  sarin  at  the 
HF  /  6-311++G(2d.2p)  level  of  theory  and  basis  set  was  approximately  43  days  of 
CPU  time  per  processor  using  the  aforementioned  hardware.  Computation  time  can 
be  reduced  by  using  several  processors  in  parallel.  In  the  case  of  “non-HCAO”  calcu¬ 
lations,  only  13  x  10  =  130  ab  initio  calculations  need  be  carried  out  per  conformer 
as  pairwise  interactions  are  neglected. 
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Figure  1:  The  vapour  phase  structures  of  the  three  lowest  energy *  con  formers  of  sarin ( 
(CH3)2CH0P(0)(F)CH^  calculated  at  HF  /  6-311+f-G(2d.2p).  The  S  enantiomer 
of  sarin  is  shown .  The  lowest  energy  con  former  is  denoted  L  the  next  lowest  en¬ 
ergy  confonner  is  denoted  If  Conformer  HI  lies  significantly  higher  in  energy  than 
conformers  I  and  II  (see  Table  1). 

3  Results  and  discussion 

3.1  Conformer  distributions  and  geometries 

The  geometries  of  the  three  lowest  energy  vapour  phase  conformers  of  sarin  were 
determined  at  HP  /  G“311++G(2d,2p)  and  were  found  to  be  in  agreement  with  the 
previous  work  of  Kaczmarek  et  ai  [44; .  The  absences  of  imaginary  frequencies  in 
ab  initio  frequency  calculations  have  confirmed  that  these  three  structures  represent 
minima  on  the  PES.  These  conformers  arc  presented  in  Figure  1.  The  lowest  energy 
conformer  of  sarin  is  denoted  sarin  I,  the  second  lowest  energy  conformer  is  denoted 
sarin  IL  and  the  subsequent  conformer  is  denoted  sarin  III. 

The  methyl  group  attached  to  the  phosphorous  atom  in  Figure  I  is  referred  to  as 
It  contains  three  CH  oscillators  which  are  referred  to  as  CH*,  CIE,  and  CH3. 
The  oscillator  which  lies  approximately  in  the  plane  of  the  page  is  CHp  That  oscillator 
which  points  into  the  plane  of  the  page  is  CH^  while  the  oscillator  which  points  out 
of  the  plane  is  CH3.  The  isopropyl  methyl  group  which  is  approximately  oriented 
towards  the  top  left-hand  side  of  the  page  is  referred  to  as  CH3 j  and  comprises  CH4* 
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CHg,  and  CH6.  The  oscillator  CH4  is  that  which  is  oriented  towards  the  upper  right- 
hand  side  of  the  page.  The  CHS  oscillator  is  oriented  out  of  the  plane  of  the  page 
while  oscillator  CH6  is  oriented  into  the  plane  of  the  page.  The  remaining  isopropyl 
methyl  group,  CH3,7,  consists  of  CH7,  CHg,  and  CHg.  The  CH;  oscillator  is  that 
which  is  oriented  downwards  and  into  the  plane  of  the  page.  The  orientation  of  CHS 
is  out  of  the  plane  of  the  page  and  towards  the  left  while  oscillator  CHg  is  oriented  out 
of  the  page  plane  and  towards  the  right-hand  side.  The  isolated  CH  oscillator  which 
contains  the  secondary  carbon  in  the  isopropyl  group  is  referred  to  as  CHi0.  There 
are  two  oxygen  atoms  in  sarin:  Oj  denotes  that  which  is  bound  to  both  carbon  and 
phosphorous  while  02  denotes  that  atom  which  is  bound  to  phosphorous  only.  The 
single  geometric  parameter  which  most  clearly  distinguishes  the  sarin  conformers  is 
the  dihedral  angle  (through  the  molecular  skeleton)  between  the  phosphorous  atom 
and  the  H  atom  in  CH|0l  denoted  Hio-  The  values  of  this  characteristic  angle,  ZH[0- 
C-Oi-P,  are  presented  for  the  S  enantiomers  of  the  three  sarin  conformers  in  Table  1. 


Table  1:  Ab  initio-calculated  dihedral  angles  and  relative  en¬ 
ergies  of  the  three  lowest  energy  vapour  phase  conformers  of 
sarin  calculated  at  HF  /  6-311++G(2d,2p) 


CONFORMER 

PARAMETER 

I 

11 

III 

DIHEDRAL  ANGLE 

ZH10-C-Oi-Pa 

-28° 

18° 

168° 

RELATIVE  ENERGIES’5 

A  E 

0  cm-1 

46.6  cm-1 

675  cm-1 

AEzpve 

0  cm"1 

29.7  cm"1 

719  cm-1 

8  This  parameter  pertains  to  the  S  enantiomer  of  sarin. 
b  Both  ZPVE  corrected  energies,  AEzpve,  and  non-ZPVE 
corrected  energies,  A E.  are  presented. 


Previously,  van  der  Waals  radii  have  been  used,  in  conjunction  with  ab  initio  calcu¬ 
lated  geometries,  as  crude  metrics  to  quantify  the  degree  of  steric  crowding  within  a 
molecule  [33].  The  internuclear  separations  between  hydrogen  H,  and  an  electronega¬ 
tive  atom  Y  (in  sarin  V-F  or  O)  to  which  it  is  not  bound,  can  be  denoted  dH,  •  Y. 
At  internuclear  separations  dH,  ■  ■  •  Y  of  approximately  2.6  A  or  less  the  separation 
between  atoms  H  and  Y  is  sufficiently  small  that  nonbonded  interactions  through 
space  are  likely  to  perturb  the  vibrational  overtone  spectrum  of  a  molecule.  Such 
through-space  interactions  are  likely  to  occur  over  internuclear  separations  larger  than 
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approximately  2,6  A  in  cases  where  large  amplitude  stretching  motions  (which  accom¬ 
pany  vibrational  overtone  transitions)  are  accompanied  by  methyl  libration  [45,46]. 
One  consequence  of  through-space  coupling  is  that  the  initially  prepared  vibrational 
stretching  overtone  state  rapidly  de-excites  via  intramolecular  vibrational  energy  re¬ 
distribution  (IVR)  [12,47],  The  spectral  manifestation  of  IVR  is  peak  broadening 
caused  by  lifetime  uncertainty  which  results  in  a  Loren tzian  profile.  Based  upon  the 
similarity  between  the  mternuciear  separations  dH,  ■  ■  ■  Y  in  sarin  [48]  and  those  in  the 
tert-butyl  halides  [33] ,  reasonable  full- width  at  half-maximum  (FWHM)  estimates  for 
the  Lorentziari  profiles  of  the  vibrational  overtone  transitions  of  sarin  are  40  cm-1  at 
Auch  =  3,  60  cm"1  at  Auch  —  4,  and  70  cm"1  at  A^ch  =  5,6*  The  HCAO  vibra¬ 
tional  overtone  spectral  profiles  of  sarin  have  been  simulated  using  slightly  smaller 
FWHM  values  than  those  indicated  above  so  as  not  to  obfuscate  detail  unnecessarily* 

The  zero-point  vibrational  energy  (ZPVE)  corrected  energies  (at  298*15  K  and  LOO  atm 
pressure)  relative  to  sarin  I,  AEzpyg,  have  been  calculated  at  HF  /  6-311-b+G(2d,2p) 
and  are  presented  in  Table  1,  along  with  the  non-ZPVE  corrected  relative  energies, 
A E.  The  AEzpve  value  for  sarin  II  was  determined  to  be  29*7  cm'1,  which  is  in  good 
agreement  with  that  reported  by  Kaczmarek  et  ai  [44],  who  determined  the  energy 
difference  to  be  31*48  cm-1  at  MP2  /6-311++G(d,p).  The  agreement  between  the 
sarin  III  AEzpve  in  Table  1  (i.  e.,  719  cm"1)*  and  that  reported  by  Kaczmarek  et 
ai  [44]  ( i .  e,,  454*68  cm”1)  is  unexpectedly  poor*  However,  Right  Walker  et  ai  [49] 
reported  the  sarin  III  AEzpve  to  be  604.5  cm-1  at  HF  /  6-311G{d,p)  and  500.8  cm"1 
at  MP2  /  6-31lG(d,p).  The  inclusion  of  some  electron  correlation  via  MP2  theory  is 
seen  to  have  had  a  pronounced  effect  on  the  relative  energy  of  the  sarin  III  confonner* 
The  relative  energies  of  the  sarin  conformers  in  Table  1  have  been  used  to  calculate 
their  relative  populations  at  298  K  in  order  to  scale  their  relative  contributions  to 
the  room  temperature  overtone  spectrum  of  sarin.  Based  upon  their  small  difference 
in  energy,  a  nearly  racemic  mixture  of  sarin  I  and  II  will  be  present  at  room  temper¬ 
ature.  The  high  relative  energy  of  sarin  III  is  such  that  this  confonner  will  have  a 
negligible  influence  and  its  spectral  contributions  are  neglected  in  this  report. 

3.2  Local  mode  and  HCAO  parameters 

Low  barriers  to  internal  methyl  rotation  (i.  e.,  up  to  a  few  tens  of  wavenumbers)  can 
allow  the  possibility  of  hundreds  of  transitions  which  contribute  to  the  observed  spec¬ 
tral  profile  [50-52],  greatly  complicating  the  observed  overtone  spectrum.  Barriers 
to  internal  methyl  rotation  of  ca.  >  450  cm-1  will  effectively  “lock”  a  methyl  group 
into  a  position,  breaking  the  coupling  between  methyl  torsion  and  CH  stretching  mo¬ 
tions  so  that  free  methyl  rotation  will  not  complicate  observed  overtone  spectra  (53). 
The  barriers  to  internal  methyl  rotation  for  conformers  I  and  II  of  sarin  have  been 
calculated  at  HF  /  6-3ll-H-G(2d,2p)  and  are  presented  in  Table  2.  In  sarin  I  the 
barriers  to  internal  methyl  torsion  were  730  crn-’  for  CH3o ,  1310  cm'1  for  CH3  .j, 


10 


DRDC  Suffield  TR  2006-220 


and  1330  cm-1  for  CH3i7  while  those  barriers  in  sarin  II  were  690  cm-1,  1340  cm^1, 
and  1370  crn-1,  respectively.  Visualization  of  the  ab  initio  output  indicates  that  no 
major  skeletal  motions  have  accompanied  methyl  torsion.  These  barriers  are  only 
approximate:  they  have  not  been  ZPVE  corrected  and  no  scaling  factors  have  been 
applied.  Nevertheless,  these  barriers  are  sufficiently  high  that  observed  CH  stretching 
overtone  spectra  will  not  he  complicated  by  methyl  torsion  structure.  The  barrier 
to  internal  rotation  of  the  CH3.*  methyl  group  has  been  previously  estimated  to  be 
677.0(4)  cm-1  [49]  from  microwave  studies  at  2  K  where  only  the  sarin  I  conformer 
should  exist  in  appreciable  quantities. 


Table  2:  Ah  initio  barriers  to  interna]  methyl  rota¬ 
tion  (in  cm~l)  of  the  two  lowest  energy  con  formers 
of  sarin  obtained  at  HF  /  6-311-h-hG(2dt2p) 


METHYL  GROUP  BARRIER 

CONFORMER 

ch3i0  CH3i/j 

ch3(. 

I 

730  1310 

1330 

II 

690  1340 

1370 

a  Barriers  have  not  been  ZPVE  corrected.  No 
scaling  factors  have  been  applied  to  the  values 
presented  in  this  table. 


The  overtone  transition  oscillator  strengths  (equation  12)  were  calculated  using  the 
transition  dipole  moments  (equation  13)  and  Morse  wavef unctions,  which  are  defined 
by  the  vibrational  quantum  number,  v)  and  the  local  mode  frequency  and  anhar- 
monicity,  uj  and  ujz  [12].  The  ab  initio  PESs  (along  the  ten  CH  stretch  coordinates 
in  sarin)  were  calculated  at  HF  /  6-311±+G(2d,2p)  in  CH  displacement  coordinate 
q  €  [-0.20,0.20]  A  in  steps  of  0,05  A,  The  vapour  phase  parameters  u  and  u )z  of 
sarin  I  and  II  were  determined  from  nonlinear  fits  of  the  ab  initio  PESs  to  equation  1 
and  subsequently  scaled  by  =  0.9429  ±0.0015  for  uj  and  sfcx  =  0.921  ±0.015  for 
u)x  [33].  These  values  of  u)  and  Cjz  are  presented  for  the  ten  CH  oscillators  of  sarin  I 
and  II  in  Table  3. 
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Table  3:  Ab  initio  calculated,  scaled ,  vapour  phase  local  inode 
parameters”  Zj  and  u )x  (in  cm~[)  of  the  two  lowest  energy  con- 
formers  of  sarin  obtained  a£  HF  /  6-311+-hG(2d,2p) 


OSCILLATOR 

CONFORMER 

I 

II 

UJ 

LUX 

Zj 

LUX 

CH, 

3095  ±  5 

58.5  ±  1.0 

3098  ±  5 

58.3  ±  1.0 

ch2 

3102  ±  5 

58.4  ±  1.0 

3103  ±  5 

58.4  ±  1.0 

ch3 

3101  ±5 

58.5  ±  1.0 

3102  ±5 

58.3  ±  1.0 

ch4 

3079  ±  5 

59.1  ±  1.0 

3077  ±  5 

59.4  ±  1.0 

CHS 

3066  ±  5 

59.3  ±  1.0 

3067  ±  5 

59.2  ±  1.0 

ch6 

3066  db  5 

59.5  ±  1.0 

3066  ±  5 

59.5  ±  1.0 

ch7 

3065  ±  5 

59.5  ±  1.0 

3061  ±  5 

59.5  ±  1.0 

ch8 

3068  ±  5 

59.2  ±  1.0 

3064  ±  5 

59.3  ±  1.0 

ch9 

3076  ±  5 

59.5  ±  1.0 

3089  ±  5 

58.9  ±  1.0 

CH10 

3094  ±  5 

60.0  ±  1.0 

3075  ±  5 

60.7  ±  1.0 

a  Local  mode  parameters  have  been  calculated  from  a  nonlin¬ 
ear  fit  of  ab  initio  data  to  equation  1  and  scaled  using  the 
appropriate  scaling  factors  from  Reference  [33]. 


The  local  mode  coupling  parameters  in  equation  7  were  calculated  in  a  manner  idem 
tical  to  that  of  Kjaergaard  et  ai  [35]  and  are  presented  in  Table  4,  The  parameter 
was  obtained  from  an  expansion  series  fit  of  the  CH  stretching  PES  [32,35]  calculated 
ab  initio  in  the  CH  displacement  range  q  6  [—0,30,0.30]  A  in  steps  of  0.05  A.  The 
parameter  7^  was  calculated  using  equation  10  and  the  ab  initio  molecular  geometry. 
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Table  4:  Ab  initio  calculated  local  mode  coupling  parameter*? 
of  the  two  lowest  energy  conformers  of  sarin  obtained  at 
HF  /  6-311++G(2d.2p) 


OSCILLATORS 

CONFORMER 

I 

II 

7ij 

% 

filj 

% 

CHj.  CHZ 

0.00317 

0.0132 

30.9 

0.00314 

0.0132 

31.2 

CHi,  CH3 

0.00339 

0.0127 

28.8 

0.00341 

0.0127 

28.8 

ch2.  ch3 

0.00322 

0.0131 

30.6 

0.00328 

0.0130 

30.2 

ch4.  ch5 

0.00474 

0.0123 

23.2 

0.00474 

0.0121 

22.6 

ch4.  ch6 

0.00481 

0.0123 

23.1 

0.00474 

0.0123 

23.3 

ch5.  ch« 

0.00466 

0.0123 

23.5 

0.00459 

0.0124 

23.8 

ch7.  ch8 

0.00472 

0.0123 

23.3 

0.00466 

0.0123 

23.5 

ch7.  ch9 

0,00456 

0.0126 

24.7 

0.00486 

0.0124 

23.1 

CHs,  CHs 

0.00454 

0.0123 

23.8 

0.00476 

0.0122 

23.1 

a  Local  mode  coupling  parameters  were  calculated  in  a  manner  iden¬ 
tical  to  that  of  Kjaergaard  et  al.  [35]. 


3.3  Simulated  spectra 

Electromagnetic  fields  select  transitions  into  states  where  all  vibrational  energy  is 
localized  into  one  of  a  set  of  equivalent  XH  oscillators.  In  the  case  of  a  methyl  group 
with  local  Cj  symmetry  {e.  <j..  sarin),  the  dominant  overtone  transitions  will  be  out  of 
the  ground  state  intestates  |i*i}|0}|0}7  10)1^)10),  or  |0)|0)|ujt).  Such  states  are  referred 
to  as  pure  local  mode  states  since  all  vibrational  excitation  is  localized  in  an  individ¬ 
ual  oscillator.  Local  mode/local  mode  combination  states  of  the  form  \vt  —  l)|l}|0), 
\vt  “  l)|0)|l),  etc.7  can  couple  into  the  pure  local  mode  states  and  steal  intensity. 
At  a  given  vibrational  overt  one ,  transitions  to  pure  local  mode  states  occur  at  lower 
frequencies  than  do  transitions  to  local  mode/local  mode  combination  states  owing 
to  oscillator  anharraonicity.  The  separation  in  the  frequencies  of  transition  to  pure 
local  mode  states  vs.  transitions  to  local  mode/ local  mode  combination  states  grows 
with  increasing  vibrational  excitation,  while  the  intensities  of  transitions  to  local 
mode/local  mode  combination  states  decrease  (as  compared  to  transitions  to  pure  lo¬ 
cal  mode  states)  with  increasing  vibrational  excitation.  The  HC AO-calculated  vapour 
phase  vibrational  overtone  transition  frequencies  and  intensities  in  the  Avch  —3  —  6 
regions  are  presented  in  Tables  B.l  and  B.2  for  a  Boltzmann-distributed  population  of 
sarin  1  and  II  at  298  K.  Individual  HC  AO-calculated  transitions  in  the  Avcu  —  3  —  6 
regions  appear  as  vertical  lines  in,  respectively.  Figures  2  through  5.  The  line  positions 
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denote  transition  frequencies  and  their  heights  indicate  relative  individual  intensities 
Each  individual  transition  has  been  assigned  a  Lorentzian  profile  with  a  FWHM  ap¬ 
propriate  to  the  degree  of  vibrational  excitation  (L  30  cm-1  at  Ahch  —  3.  40  cm" 1 
at  Avcu  —  4.  and  50  cm-1  at  Avch  —  5,6}  and  the  individual  Lorentzian  profiles 
have  been  summed  to  yield  HC AO-simulated  spectra.  It  is  noteworthy  that  while  the 
dominant  transitions  in  the  Ai^ch  =  3  —  6  regions  are  transitions  into  pure  local  mode 
states,  combination  hands  make  significant  contributions  to  the  observed  vibrational 
spectra,  especially  at  lower  overtones  (see  below). 


Figure  2:  The  vapour  phase  A*  vh  =  3  HCAO  (black)  and  "non-HCAQ  (rod)  cal¬ 
culated  vibrational  overtone  spectral  region  of  an  equilibrium  distribution  of  sarin  I 
and  II  eonformers  at  298  K.  The  maximum  (summed)  HCAO  absorbance  feature  in 
this  region  has  an  oscillator  strength  of  2.47  x  10 “9  at  cn.  8470  rm"1.  The  %on- 
HCAO "  Irat  e  has  been  scaled  so  that  its  maximum  ahsurhaino  feature  has  the  same 
amplitude  as  the  corresponding  feature  in  the  HCAO  trace. 


Computationally  less  expensive  “non- HCAO1’  transitions  and  spectra  have  been  sim¬ 
ulated  in  the  Ai'ch  —3  —  6  regions  and  are  presented  in,  respectively.  Figures  2 
through  5.  As  noted  previously,  the  intensity  data  from  ‘'non-HCAO”  calculations 
are  notoriously  unreliable  the  intensities  of  the  “non- HCAO  transitions  have  been 
scaled  so  that  the  global  maxima  of  the  summed  “non-HCACT  spectral  profiles  are 
identical  to  those  of  the  H C A O-calc u  1  ated  profiles.  Small  variations  in  the  heights  of 
individual  “non- HCAO"  transition  lint's  reflect  the  relative  abundances  of  the  sarin  1 
and  II  eonformers  at  298  K. 


14 


DRDC  Suffield  TR  2006-220 


Figure  3:  The  vapour  phase  AtfcH  =  4  HCAO  (black)  anti  “ non-HCAO ”  (red)  cal¬ 
culated  vibrational  overtone  spectral  region  of  an  equilibrium  distribution  of  sarin  I 
and  II  con  formers  at  298  K,  The  maximum  (summed)  HCAO  absorbance  feature  in 
this  region  has  an  osciiiator  strength  of  2.56  x  10“ 10  at  ca.  11060  cm'1.  The  ‘'non- 
HCAO"  trace  has  been  scaled  so  that  its  maximum  absorbance  feature  has  the  same 
amplitude  as  the  corresponding  feature  in  the  HCAO  trace. 


The  vapour  phase  HCAO  /<*«.  values  of  the  individual  Boltzmann  weighted  transitions 
within  the  v  -  3  manifold  are  shown  in  Figure  2  as  black  vertical  lines.  The  summed 
HCAO  A'Och  =  3  profile  is  shown  in  black  and  has  a  maximum  /osc  of  2.47  x  10-9 
at  ca.  8470  cm”1.  The  Boltzmann  weighted  "non-HCAO”  transitions  are  denoted  by 
red  vertical  lines.  The  “non-HCAO”  profile  is  erroneously  blue-shifted  relative  to  the 
HCAO  profile  as  the  former  neglects  interoscillator  coupling.  However,  interoscillator 
coupling  decreases  with  increasing  vibrational  excitation  [17]  so  that  the  disagreement 
between  HCAO  and  “non-HCAO”  transition  frequencies  will  decrease  with  increasing 
excitation. 

Local  inode/local  mode  combination  band  structure  is  absent  in  the  “non-HCAO” 
profile  as  the  “non-HCAO”  approach  neglects  interoscillator  coupling.  The  presence 
of  significant  local  mode/ local  mode  combination  band  structure  above  ca.  8625  cm'1 
in  the  HCAO  simulation  indicates  that  the  accurate  calculation  of  the  vibrational 
overtone  spectrum  of  sarin  at  Al'ch  =  3  requires  the  treatment  of  inter-oscillator 
coupling.  Agreement  between  the  “non-HCAO"  and  HCAO  calculated  spectral  re¬ 
gions  is  poor,  indicating  that  the  less  expensive  “non-HCAO”  approach  is  unsuitable 
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Figure  4:  The  va pour  phase  Auqh  —  5  HCAO  (black)  and  “ non-HCAO "  (red)  cal¬ 
culated  vibrational  overtone  spectral  region  of  an  equilibrium  distribution  of  sarin  ! 
and  1 1  conformers  at  298  K.  The  maximum  (summed)  HCAO  absorbance  feature  in 
this  region  has  an  oscillator  strength  of  2.62  x  IQ-11  at  ca .  13540  cnr1.  The  “ non - 
HCAO"  trace  has  been  scaled  so  that  its  maximum  absorbance  feature  has  the  same 
amplitude  as  the  corresponding  feature  in  the  HCAO  trace . 


at  Ai^ch  —  3. 

The  simulated  HCAO  spectral  profile  at  Aren  —  3  (and,  in  fact*  at  all  overtones 
presented  in  this  report)  is  deceptively  simple  as  it  arises  from  numerous  unresolved 
CH  overtone  transitions*  Individual  vibrational  overtone  transitions  of  CH  bonds  can 
be  resolved  from  one  another  in  room-temperature  vibrational  overtone  spectra  if  the 
bond  lengths  differ  by  as  little  as  1  mA  [7,54).  However,  in  sarin  I  the  CH  bond 
lengths,  calculated  at  HF  /  6-311++G(2d,2p),  range  from  1,079  A  to  1.083  A.  with 
the  difference  between  the  most  similar  bond  lengths2  being  0.01  mA.  Similarly,  in 
sarin  II,  the  bond  lengths  range  from  1.080  A  to  1,083  A  with  a  difference  between 
the  most  similar  bond  length  values  of  0.02  mA.  A  mixture  of  sarin  1  and  II  results 
in  a  range  of  CH  bond  lengths  where  the  difference  between  the  most  similar  bond 
length  values  is  0.0004  m A.  The  CH  vibrational  overtone  transitions  in  sarin  will  not 

2 Ab  mi/ie-caleulated  bond  length  differences  on  the  order  of  0.01  m A  or  less  are  not  usually 
considered  to  be  physically  significant  as  limitations  in  theory  level  and  basis  set  are  expected  lo 
contribute  errors  which  are  similar  in  magnitude  for  all  but.  the  highest  level  calculations  carried  out 
on  relatively  small  molecules. 
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Figure  5:  The  vapour  phase  Auch  =  6  HCAO  (block)  and  “non-HCAO"  (red)  cal¬ 
culated  vibrational  overtone  spectral  region  of  an  equilibrium  distribution  of  sarin  I 
and  11  conformers  at  298  K.  The  maximum  (summed)  HCAO  absorbance  feature  in 
this  region  lias  ari  oscillator  strength  of  3.42  x  10“ 12  at  ca.  15890  cm~l.  The  “ non - 
HCAO"  trace  has  been  scaled  so  that  its  maximum  absorbance  feature  has  the  same 
amplitude  as  the  corresponding  feature  in  the  HCAO  trace. 


be  experimentally  resolvable  at  room  temperature. 

The  HCAO  and  “non-HCAO”  Ai'ch  —  4  calculated  vapour  phase  spectral  regions  are 
presented  in  Figure  3  as  individual  transitions  and  as  summed  profiles.  The  maximum 
HCAO  /(**.  is  2.56  x  10-10  and  occurs  at  ca.  11060  cm-1.  Local  mode/local  mode 
combination  bands  are  visible  above  11350  cm-1.  However,  these  combination  bands 
are  less  prominent  than  at  Auch  —  3  and  are  increasingly  becoming  separated  from 
the  pure  local  mode  bands.  The  numerous  unresolved  transitions  which  form  thp 
summed  HCAO  profile  reinforce  the  necessity  to  carry  out  HCAO  calculations  if  one 
is  to  assign  and  interpret  observed  overtone  transitions.  The  disparity  between  the 
HCAO  profile  (which  required  1534  ab  initio  calculations  per  conformer)  and  "non- 
HCAO”  profile  (which  required  130  ab  initio  calculations  per  conformer)  indicates 
that  an  HCAO  approach  is  required  to  simulate  the  Auch  =  4  spectral  region  of 
sarin. 

The  agreement  between  the  HCAO  and  “non-HCAO”  Auch  —  5  simulated  vapour 
phase  spectra,  presented  in  Figure  4,  is  passable.  Unfortunately,  at  this  degree  of 
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vibrational  excitation  (and  higher)  the  inherent  intensities  of  the  vibrational  overtone 
transitions  are  very  weak  and,  in  general,  are  no  longer  offset  by  increases  in  portable 
laser  power  and  detector  efficiency  in  the  visible  region  of  the  spectrum.  The  Ai>ch  = 
5  HCAO  maximum  summed  flKjF  is  2.62  x  10-11  at  ca.  13540  cm-1.  Local  mode/local 
mode  combination  bauds,  which  occur  above  13900  cm-1,  are  continuing  to  decrease 
in  significance  with  increasing  vibrational  excitation. 

The  HCAO  and  “non-HCAO"  simulated  vapour  phase  spectra  of  the  Ai’ch  =  0 
region  are  presented  in  Figure  5.  The  maximum  summed  At.’cn  =  0  HCAO  is 
3.42  x  10“ 12  at  ca.  15890  cm-1.  Agreement  between  the  “non-HCAO”  and  HCAO 
simulated  spectra  are  adequate  for  most  detection  (but  not  identification)  purposes. 
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4  Conclusions 


The  three  lowest  energy  vapour  phase  conformers  of  sarin  have  been  identified.  The 
two  lowest  energy  conformers  were  determined  to  contribute  significantly  to  the  spec¬ 
tral  profile  of  a  sample  of  sarin  at  289  K.  The  barriers  to  internal  methyl  rotation 
in  sarin  were  found  to  be  sufficiently  high  that  the  coupling  of  methyl  torsion  to  CH 
stretching  will  not  complicate  the  observed  overtone  spectral  regions. 

The  vapour  phase  local  mode  parameters  lD  and  u)x  have  been  calculated  for  each  of 
the  ten  CH  oscillators  in  the  two  spectrally  significant  sarin  conformers  ab  initio  at 
the  HF  /  6-311  +  +G(2d.2p)  level  of  theory  and  basis  set,  as  have  the  inter-oscillator 
coupling  parameters.  In  conjunction  with  dipole  moment  functions  derived  from 
ab  initio  calculations,  those  parameters  above  were  used  to  perform  harmonically 
coupled  anharmonic  oscillator  (HCAO)  calculations.  The  HCAO  calculations  were 
in  turn  used  to  simulate  the  vapour  phase  vibrational  overtone  spectral  regions  in 
a  room-temperature  sample  of  sarin.  It  has  been  shown  to  be  feasible  to  predict 
(without  recourse  to  the  experimental  spectra)  the  absorption  spectra  of  species  which 
are  difficult  to  synthesize,  handle,  or  otherwise  acquire. 

The  inclusion  of  inter-oscillator  coupling  (uia  the  HCAO  approach)  has  been  shown 
to  be  necessary  to  predict  the  lower  vibrational  overtone  regions  (t.  e.,  first  to  third 
overtones)  as  a  simpler  “non-HCAO”  approach  (which  does  not  allow  pairwise  har¬ 
monic  coupling  among  adjacent  oscillators)  failed  to  accurately  reproduce  the  HC AO- 
simulated  pure  local  mode  spectral  regions. 
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5  Future  work  and  recommendations 


Further  validation  of  the  approach  to  spectral  simulation  used  in  this  report  should 
be  undertaken  by  acquiring  the  experimental  overtone  spectral  regions  of  sarin  in  the 
vapour  and  liquid  phases.  A  comparison  of  the  theoretical  and  experimental  vapour 
phase  spectra  will  indicate  whether  Fermi  resonances3  are  likely  to  complicate  the 
overtone  spectra  of  organophosphate  nerve  agent  homologues  (i.  e.,  R0P(0)(X)R'), 
making  the  task  of  simulating  their  vibrational  overtone  spectra  more  challenging. 
Future  experimental  and  theoretical  investigations  into  the  vapour  phase  overtone 
spectral  regions  of  other  nerve  agents  should  be  undertaken  in  order  that  possibility 
of  deriving  local  mode  scaling  factors  specific  to  R0P(0)(X)R'  homologues  can  be 
investigated.  Finally,  the  experimental  overtone  spectra  of  nerve  agents  should  be 
acquired  in  the  liquid  phase  and  isolated  in  matrices  so  that  the  effects  of  solvation 
on  experimental  local  mode  parameters  and  overtone  transition  intensities  can  be 
compared  with  results  from  HCAO  calculations  based  upon  ab  initio  solvation  models 
in  order  to  investigate  our  capability  to  predict  the  overtone  spectra  of  CWAs  in 
condensed  phases. 

The  a b  initio  calculations  which  were  used  in  the  research  presented  in  this  report 
were  carried  out  on  a  publicly-funded  academic  computer  network.  It  would  be  ad¬ 
vantageous  for  Defence  Research  and  Development  -  Suffield  to  acquire  the  necessary 
hardware  and  software  in  order  to  be  able  to  carry  out  such  research  in  a  secure 
environment. 


3Energy  flow  between  a  pure  local  mode  state  and  near- resonant  combination  states  has  been 
observed  in  some  molecules.  Such  an  energy  flow  causes  I  lie  wavefn  net  ions  of  the  corresponding 
states  to  mix.  This  mixing  can  be  observed  as  a  splitting  in  the  frequency  domain  \A1\  which  is 
often  referred  to  as  Fermi  resonance  [55 j. 
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Annex  A:  Ab  initio  input  matrices  and 
geometries 

The  Gaussian  03  input  z- matrices  used  to  calculate  the  harmonic  frequencies  of  the 
three  lowest  energy  conformers  of  sarin  are  presented  below.  The  parameters  which 
begin  with  an  V  denote  bond  lengths  (in  A),  the  parameters  which  begin  with  an  “a” 
denote  bond  angles  (in  degrees),  and  the  parameters  which  begin  with  an  “d”  denote 
dihedral  angles  (in  degrees).  These  parameters  are  taken  from  ab  initio  geometry 
optimizations  calculated  at  HP  /  6-311++G(2d,2p)  calculated  with  SCF=Tight  and 
Opt=Tight  options  specified.  Geometry  optimizations  were  run  using  all  Gaussian 
defaults  save  that  Gaussian  overlay  option  IOP  (3/32=2)  was  used  to  prevent  the 
possible  reduction  of  the  expansion  set. 


A.t  Sarin  conformer  I 

'/,chk=If  req .  chk 

#  HF/6-311++G(2d ,2p)  Freq  SCF=Tight  IOP (3/32=2) 
Freq  (HF/6-311++G(2d,2p))  of  conformer  I  of  GB 


0  1 
P 

01  P  rOlP 

C2  01  rC20  P  aCOP 

C3  C2  rC3C  01  aC3C0  P  dC3P 

C4  C2  rC4C  01  aC4C0  C3  dC4C3 

H10  C2  rC4H  01  aHCO  C3  dHC3 

H4  C3  rH4C  C2  aH4C  01  dH40 

95  C3  rH5C  C2  aH5C  H4  dH5H4 

96  C3  rH6C  C2  aH6C  94  dH6H4 

H7  C4  rH7C  C2  aH7C  01  dH70 

H8  C4  rH8C  C2  aH8C  H7  dH8H7 

H9  C4  rH9C  C2  aH9C  H7  dH9H7 

F  P  rFP  01  aFOl  C2  dFC2 

02  P  r02P  01  aD201  F  d02F 
Cl  P  rClP  01  aClOl  F  dCIF 
HI  Cl  rHIC  P  aHIP  F  dHlF 
H2  Cl  rH2C  P  aH2P  HI  dH2Hl 
H3  Cl  rH3C  P  aH3P  HI  dH3Hl 

rOlP  =  1.55219807 
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rC20  =  1.44561224 
rC3C  =  1.51545258 
rC4C  =  1.5132104 
rC4H  =  1.07882276 
rClP  =  1.78014048 
rFP  =  1.54526258 
r02P  =  1.43778121 
rHIC  =  1.08069984 
rH2C  =  1.08005479 
rH3C  =  1.08023 
rR4C  =  1.0815011 
rH5C  =  1.08288668 
rH6C  =  1.082665 
rH7C  =  1.08272136 
rH8C  =  1.08270156 
rH9C  =  1.0814916 
aCOP  =  123.8916557 
aC3C0  =  109.14458558 
aC4C0  =  106.58583333 
aHCO  =  107.3419504 
aClOl  =  104.44852 
aFOl  =  102.69551244 
a0201  =  116.30956285 
aHIP  =  108.81848874 
aH2P  =  110.35815568 
aH3P  =  108.96179048 
aH4C  =  111.07003332 
aH5C  =  109.75615367 
aH6C  =  110.60899134 
aH7C  =  110.67406835 
aHSC  =  109.8723145 
aH9C  =  110.62820183 
dC3P  =  91.12722011 
dC4C3  =  122.85968156 
dHC3  =  -119.29976504 
dFC2  =  -96.07838913 
d02F  =  122.49979869 
dCIF  =  -106.1888713 
dH4Q  =  -61.11664106 
dH5H4  =  -119.57527707 
dH6H4  =  120.60261824 
dH70  =  -61.19614335 
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dH8H7  =  -119.86383057 
dH9H7  =  120.44445118 
dHlF  =  63.93484782 
dH2Hl  =  -120.6895776 
dH3Hl  =  118.905770 


A.2  Sarin  conformer  II 

*/,chk=IIfreq.  chk 

#  HF/ 6-31 1++G (2d , 2p)  Freq  SCF=Tight  I OP (3/32=2) 
Freq  (HF/6-311++G(2d,2p) )  of  conformer  II  of  GB 


0  1 

P 

01  P  rOlP 
C2  01  rC20  P  aCOP 
C3  C2  rC3C  01  aC3C0  P  dC3P 
C4  C2  rC4C  01  aC4C0  C3  dC4C3 
H10  C2  rC4H  01  aHCO  C3  dHC3 
H4  C3  rH4C  C2  aH4C  01  dH40 

H5  C3  rH5C  C2  aH5C  H4  dH5H4 

H6  C3  rH6C  C2  aH6C  H4  dH6H4 

H7  C4  rH7C  C2  aH7C  01  dH70 

H8  C4  rH8C  C2  aH8C  H7  dH8H7 

H9  C4  rH9C  C2  aH9C  H7  dH9H7 

F  P  rFP  01  aFOl  C2  dFC2 
02  P  r02P  01  a0201  F  d02F 
Cl  P  rClP  01  aClOl  F  dCIF 
HI  Cl  rHIC  P  aHIP  F  dHlF 
H2  Cl  rH2C  P  aH2P  HI  dH2Hl 
H3  Cl  rH3C  P  aH3P  HI  dH3Hl 

rOlP  =  1.55209928 
rC20  =  1.44687758 
rC3C  =  1.51356472 
rC4C  =  1.51379742 
rC4H  =  1.08005516 
rClP  =  1.77979872 
rFP  =  1.54609944 
r02P  =  1.43770942 
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rHIC  =  1.08071916 
rH2C  =  1.07995462 
rH3C  =  1.08023872 
rH4C  =  1.08146989 
rHSC  =  1.08282074 
rH6C  =  1.08264867 
rH7C  =  1.08302149 
rH8C  =  1.0830462 
rH9C  =  1.08064641 
aCOP  =  122.88457571 
aC3C0  =  106.9581287 
aC4C0  =  108.45230988 
aHCO  =  108.10681882 
aClOl  =  104.24119054 
aFOl  =  102.54531762 
a0201  =  116.52542919 
aHIP  =  108.75430263 
aH2P  =  110.39447289 
aH3P  =  108.96024269 
aH4C  =  110.69797645 
aH5C  =  109.89538515 
aH6C  *  110.56103064 
aH7C  =  110.45813428 
aH8C  =  109.89473145 
aH9C  =  110.44580844 
dC3P  =  135.99520005 
dC4C3  =  123.12764795 
dHC3  =  -118.26612432 
dFC2  =  -87.52650936 
d02F  =  122.14272659 
dCIF  =  -106.4215882 
dH40  =  -59.45350453 
dH5H4  =  -119.84265928 
dH6H4  =  120.35992967 
dH70  =  -60.05386025 
dH8H7  =  -119.69636064 
dH9H7  =  120.66341364 
dHlF  =  64.38551409 
dH2Hl  =  -120.55506988 
dH3Hl  =  118.83160012 
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A.3  Sarin  conformer  III 


%chk= I I I f r eq . chk 

#  HF/6-311+-*-G(2d,2p)  Freq  SCF=Tight  I  OP  (3/32=2) 

Freq  (HF/6-311++G(2d,2p))  of  conformer  III  of  GB 

0  1 
P 

01  P  rOlP 
C2  01  rC20  P  aCOP 
C3  C2  rC3C  01  aC3C0  P  dC3P 
C4  C2  rC4C  01  aC4C0  C3  dC4C3 
H10  C2  rC4H  01  aHCO  C3  dHC3 
H7  03  rH7C  C2  aH7C  01  dH70 

H8  C3  rH8C  C2  aH8C  H7  dH8H7 

H9  C3  rH9C  C2  aH9C  H7  dH9H7 

H4  04  rH4C  02  aH4C  01  dH4Q 

H5  04  rH5C  02  aHSC  H4  dH5H4 

H6  04  rH6C  02  aH6C  H4  dH6H4 

F  P  rFP  01  aFOl  02  dFC2 
02  P  r02P  01  a0201  F  d02F 
Cl  P  rClP  01  aClOl  F  dCIF 
HI  Cl  rHIC  P  aHIP  F  dHlF 
H2  Cl  rH2C  P  aH2P  HI  dH2Hl 
H3  Cl  rH3C  P  aH3P  HI  dH3Hl 

rOlP  =  1 . S5022497 
rC20  =  1.44762438 
rC3C  =  1.51575231 
rC4C  =  1.51693348 
rC4H  =  1.07881959 
rClP  =  1.78063789 
rFP  =  1.54576025 
rD2P  =  1.43810658 
rHIC  =  1.08074166 
rH2C  =  1.07991007 
rH3C  =  1.08029348 
rH7C  =  1.08213067 
rH8C  =  1.08325202 
rH9C  =  1 . 08008308 
rH4C  =  1.08150232 
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rH5C  =  1.08272806 
rH6C  =  1.08157104 
aCOP  =  127.93600428 
aC3C0  =  110.05938518 
aC4C0  =  111.50018795 
aHCO  =  102.38426845 
aClOl  =  103.99577746 
aFOl  =  102.60792751 
a0201  =  117.30057445 
aHIP  =  108.74024546 
aH2P  =  110.44738668 
aH3P  =  108.97646918 
aH7C  =  110.39311746 
aH8C  =  109.33059246 
aH9C  =  111.03540767 
aH4C  =  111.31335865 
aH5C  =  109.06713619 
aH6C  =  110.78627298 
dC3P  =  -76.26372773 
dC4C3  =  127.28586373 
dHC3  =  -116.21320077 
dFC2  =  -92.58859766 
d02F  =  122.70973745 
dCIF  =  -106.06245935 
dH70  =  -53.63459719 
dH8H7  =  -118.97432431 
dH9H7  =  120.85933206 
dH40  =  -68.05096474 
dH5H4  =  -119.4588776 
dH6H4  =  121.40441189 
dHlF  =  64.51430506 
dH2Hl  =  -120.57938587 
dH3Hl  =  118.81112983 
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Annex  B: 

Calculated 

oscillator  strenaths 

Table  B.1: 

Transition  frequencies  (in  cm  1  )>  i/v 

,—0,  and  oscillator  strengths ,  fQ&C}  of 

the  HCAO  calculated*  vapour  phase 

CH  vibrational  overtone  transitions  An  =  3  —  6 

in  sarin  I 

|CH1)|CH1)|CH3) 

tCH4)|CH5)|CH6) 

|CH7)]CH8)|CH9) 

|CH10) 

State 

0 

faze 

/oBC 

/(lie 

Uu*—  O  lose 

!3>|0)|0) 

8564 

6.99  x  10- 10 

8502 

7.35  X  10" 10 

8455 

7.06  x  10“ 10 

|2)|1)|0) 

8794 

5.74  x  10“U 

8806 

4.63  x  10“12 

8650 

2.08  x  10- 10 

|1)|2)|0) 

8901 

l  20  X  10“ 10 

8705 

8  25  x  10-11 

8768 

6.61  x  10"n 

|0>|3>|0> 

8579 

5,07  x  10“10 

3475 

9.03  x  10“10 

8467 

8.59  x  10" 10 

|a>|o>|i> 

8798 

1,73  x  10“10 

8784 

6.86  x  10“ 11 

8704 

1,12  x  10“ 10 

|o)W|i> 

8927 

7.00  x  10“13 

8776 

9.84  x  10“ 11 

8715 

3.10  x  10"10 

twm 

8S94 

1.39  x  10"11 

8698 

2.79  x  10“ 10 

8799 

4,72  x  10“n 

[0>[1)|2> 

8736 

2.68  x  10“ lu 

8652 

2.08  x  10" 10 

8823 

1,78  x  1G“U 

|0)|0)|3> 

8585 

5.31  x  10-'° 

8468 

6.81  x  10" 10 

8544 

8.79  x  10“in 

Au  =  3  total 

2.38  x  10“® 

3,07  x  IQ“9 

3.22  x  i0“® 

8497  1,10  x  10“® 

14)10)10} 

11203 

6.20  x  10“ 11 

11107 

8.36  X  10" 11 

11042 

6,77  x  10“  H 

|3>|1>|0> 

11474 

4.78  x  10“ 13 

11477 

7.54  x  10"13 

11360 

3.37  x  10“12 

|t>|3>|0> 

11587 

7.70  x  10“ 12 

11402 

1.23  x  10-11 

11400 

2.04  x  10“n 

j0)|4>[0) 

11221 

6.02  x  10“ 11 

11069 

8.70  x  10“  11 

11058 

8.62  x  10"11 

|3)  |0)|1) 

11523 

3.43  x  I0"ia 

11449 

2.58  x  IQ" 12 

11419 

2.17  x  10"12 

!0}|3)[1) 

11526 

1,23  x  10"  1 1 

1 1367 

4,60  x  10“12 

11442 

4,24  x  10“12 

|1>|0>|3> 

11584 

4.81  x  10“a3 

11406 

1.28  x  I0“u 

11469 

4,48  x  W~  V2 

[0)|1)|3) 

11618 

1.12  X  1G“13 

11445 

5.48  x  10“ 12 

11507 

223  x  10“ 12 

|0)  f0)|4) 

1 1219 

4,35  x  10“13 

11061 

6,75  x  10“u 

11165 

9,45  x  IQ“n 

A  u  —  4  total 

1,97  x  10“ 10 

2,79  x  10“ 10 

2.87  x  10~ao 

}  1086  8.94  x  10” 11 

|5)|0)|0) 

13720 

6.27  x  10“ 12 

13591 

9,34  x  10“ 12 

13508 

7.05  x  10“ 12 

|4>|1>|0> 

14137 

2,55  x  10“ 13 

14071 

8.10  x  10“ 14 

13965 

3.92  x  10“i3 

|1>|4>|0> 

l  116-1 

6.98  x  10“13 

14039 

1.63  x  10“  M 

13981 

1.29  x  lQ~i2 

|0)|5)|0) 

13741 

5.69  x  10“12 

13545 

9.20  x  10“ 12 

13531 

8.89  x  10“12 

|4)|0)|!) 

14206 

1.76  x  10“ 13 

14020 

1,13  x  10“12 

14020 

1,01  x  1G-*3 

|0)|4)jl) 

14216 

1.31  x  10“ 13 

13991 

1.14  x  10“12 

14036 

2.75  x  10- 13 

|1)|0)|4) 

14231 

2,80  x  10“ 14 

14031 

2.49  x  10" 13 

14076 

3,49  x  10-*3 

|0)|I)|4) 

14159 

7.27  x  10“ 13 

13981 

2,70  x  10“ 13 

14123 

2,85  x  1Q“13 

|0)|0)|5) 

13741 

4,57  x  10“ 13 

13533 

6,94  x  10“ 12 

13667 

9.19  x  10“12 

Au  ~  5  total 

1.87  x  10“n 

2.85  x  10" 11 

2,80  x  10“ 11 

13554  8.82  x  10" 12 

|6)|0)|0) 

16120 

8.62  x  10"13 

15956 

1,36  X  10  " 12 

15856 

1.03  x  10-12 

|5>U>|Q> 

16728 

7.65  x  10"1& 

16557 

8,08  x  10“ 15 

16436 

7.51  x  ID"14 

|1}|5)|0) 

16687 

9.59  x  10"M 

16517 

5.34  x  10“15 

16455 

9.49  x  10- 14 

|0}|6)|0) 

16144 

8.98  x  10“13 

15901 

1.24  x  10“12 

15884 

1,22  x  10“ 12 

|S}|0)|1) 

16667 

6.38  x  10“14 

16509 

9,76  x  10“14 

16489 

1,44  x  10“ 14 

(D>|5>|1> 

16745 

7.44  x  10"15 

16468 

1.24  x  10“13 

16510 

2,19  x  10“ 14 

|1)|0>|5> 

166-84 

2.85  x  10“H 

16505 

2. 38  x  10“  N 

16581 

3,21  x  10“ 14 

|0>|1>|5> 

16749 

1.04  x  10“14 

16458 

2.81  x  10“14 

16629 

2,82  x  10“ 14 

|0>|0)  |6) 

16144 

7.29  x  10“13 

15886 

9,51  x  10“13 

16050 

1.28  x  10“ 12 

iu  =  6  total 

2.72  x  10“12 

3,85  x  10“12 

3.80  x  10-*2 

15901  1.10x  10" 12 

*  An  13  point  grid  was  used  to  calculate  Lhe  dipole  moment  function  derivatives  in  equation  14, 


DRDC  Suffield  TR  2006-220 


33 


Table  8.2:  Transition  frequencies  (in  cm  1),  uv^o,  and  oscillator  strengths,  fosc,  of 
the  HCAO  calculated11  vapour  phase  CH  vibrational  overtone  transitions  Av  =  3  —  6 


in  sarin  II 


State 

|CH,)|CH2>|CH3> 

iCH4>|CH5)|CH6) 

|CHT)|CHa>|CH*) 

ICHio) 

/owe 

/q&c 

/osc 

/o»c 

|3)|0>10) 

8556 

6.62  x  10- 10 

8510 

8,55  x  IO' 10 

8465 

6,63  x  10“ 10 

12)11)10) 

8789 

5.34  x  10-“ 

8809 

8.83  x  10- 12 

6649 

2,07  x  10" 10 

l»|2)|0) 

8897 

]  .14  x  icr10 

8699 

2.93  x  10- 10 

8775 

6,02  x  10“ 31 

|0)|3)|0> 

8581 

5.42  X  10- 10 

8472 

8.99  x  10-1D 

8478 

8.03  x  10- 10 

|2)|C»|t> 

8994 

1.28  X  10-»> 

8787 

8.58  x  1G-ID 

3697 

1.33  x  10“ 10 

|0)|2)|1> 

8923 

2.24  X  10- 12 

8774 

8.06  x  IO” 11 

8703 

2.79  x  10“ 10 

|l)|0)|2) 

8888 

1.77  x  10- 11 

8707 

9.66  x  10“ 11 

8786 

3.88  x  10- 3  3 

|0>|1>|2) 

8793 

1.73  x  10' 10 

8652 

2.28  x  IO" 10 

8808 

5.78  x  IO" 33 

|0>|0>|3) 

8574 

4.97  x  10" 10 

8467 

7  02  x  10“ 10 

8498 

7,29  x  10- 10 

Av  =  3  total 

2.33  x  10 -® 

3.26  x  IO”0 

2.99  x  10‘9 

8561 

1.27  x  IQ-9 

14>|0)|0> 

11189 

7.69  X  10“  u 

11119 

8.15  x  10“31 

11057 

7  16  x  10“n 

13)11)10} 

11466 

5.78  x  10- 12 

11482 

9,80  x  10“33 

11364 

3,99  x  10“ 12 

|1>|3)|0) 

11580 

7.64  x  10- 12 

11402 

1,01  x  10-11 

11399 

2.31  x  1G-U 

|0>|4)|0) 

11216 

5,24  x  10“ 11 

11065 

S.59  x  10“n 

11074 

9  98  x  10“ u 

|3)|0>|1) 

11515 

4.41  x  10“12 

11452 

5.37  x  10- 12 

11443 

1.33  x  10  12 

|0>|3)|1) 

11518 

1.10  X  io- 11 

11367 

5,36  x  10“  12 

11449 

3.03  x  10“*2 

[1  >[0>  |3) 

11576 

2.92  x  10“ 13 

11410 

1,50  x  10"n 

11404 

5.17  x  1Q“12 

[0>|1)|3> 

11611 

2.38  x  10“ 13 

11445 

2.64  x  10“32 

11475 

4.44  x  10“*2 

|0)|t>)|4) 

11213 

5.02  x  IQ-11 

11060 

7.14  x  I0“n 

11100 

8.09  x  10“ 11 

Av  =  4  total 

2.11  x  10"10 

2.80  x  IO” 10 

2,96  x  10- 30 

11175 

9.98  x  I0“n 

|5)|0)|0> 

13701 

6.43  x  10'12 

13607 

9.17  x  IQ"12 

13528 

6.86  x  10-12 

|4>|1>|0> 

14124 

3.58  x  IO-13 

14027 

8.06  x  10“ 13 

13977 

4.75  x  10“ 13 

11>|4>|0> 

14157 

6,75  x  10- 13 

14038 

3.63  x  10" 13 

13993 

1.40  x  10“12 

|0>|5)|0) 

13736 

5.34  x  IO" 12 

13539 

9.09  x  IO" 12 

13551 

8,96  x  IO'32 

[4>[0>|1> 

14193 

2,12  x  10'13 

14031 

8.90  x  10-1J 

14028 

1.77  x  IO"33 

|0)|4)H) 

14207 

1.27  x  10“ 13 

13989 

1.11  x  10“13 

14044 

8  79  x  10“ 14 

|l>t0>t4) 

14149 

5.89  x  10- 13 

14035 

2,05  x  10-13 

14014 

2,92  x  10“ 13 

|0)|1}|4) 

14223 

3.48  x  10“ 14 

13979 

3.51  x  10“ 13 

14064 

4.80  x  10- 33 

|0)|0)|S) 

13731 

4,79  x  IO" 12 

13532 

7,07  x  10” 12 

13583 

9  03  x  10- 12 

Av  =  5  total 

1.87  x  10“ u 

2.84  x  IO"31 

2,79  x  10“11 

13669 

9.46  x  10“ 12 

16)|0>|0} 

16095 

1.13  x  10- 12 

15977 

1.19  x  10“32 

15880 

1  01  x  10“ 32 

|5>|1)|0) 

16709 

8.17  x  1Q“1S 

16573 

8.88  x  10' 35 

16452 

S.29  x  10“ 14 

11)15)10} 

16679 

8.58  x  10"  14 

16512 

6.38  x  10“ 15 

16472 

1.06  x  10“ 13 

|0)|6>!0> 

16138 

7 .84  x  10- 13 

15894 

1.27  x  10“ 12 

15909 

1  30  x  10“12 

|5i|0)|l) 

16647 

8.07  x  10“u 

16525 

9.57  x  1Q“14 

16501 

1.26  x  10“15 

10)15)10 

16739 

1.82  x  IO" 15 

16463 

1,30  x  10“  13 

16524 

8.22  x  10" 35 

|1)10)]5) 

16673 

2.18  x  IO"14 

16506 

1.95  x  10“ 14 

16502 

4,62  x  HI-14 

|0)|1)[5) 

16736 

1.75  x  10“M 

16457 

3  GO  x  10“ 14 

16549 

4.96  x  10“  N 

10)10)  |6) 

16132 

7.50  x  ID”13 

15885 

9  98  x  10“13 

15946 

1.34  x  10' 32 

Av  =  6  total 

2.89  x  10“ 12 

3.76  x  10“ 12 

3.95  x  10' 12 

16043 

1.14  x  10“12 

n  An  11  point  grid  was  used  to  calculate  the  dipole  moment  function  derivatives  in  equation  14 
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In  sarin  (isopropyl  melhylphosphonoftuoridate)  there  are  ten  nonequivalent  CH  oscillators.  Ab 
initio  calculations  at  the  HF/  6-31  U+G(2d,2p)  level  have  been  used  to  determine  the  vapour 
phase  local  mode  parameters,  w  and  mx,  for  each  oscillator  in  the  two  spectrally  significant 
conformers  of  sarin,  as  well  as  inter-oscillator  coupling  parameters.  These  above  parameters, 
in  conjunction  with  dipole  moment  functions  derived  from  ab  initio  calculations,  were  used  to 
perform  harmonically  coupled  anharmonic  oscillator  (HCAO)  calculations,  thereby  enabling  the 
simulation  of  vibrational  overtone  spectral  regions  in  a  room-temperature  sample  of  sarin  It  was 
determined  that  the  computafionally-intensive  HCAO  approach  is  necessary  to  predict  the  lower 
vibrational  overtone  regions  (re.,  first  to  third  overtones)  as  a  simpler  "non-HCAO"  approach 
(which  does  not  allow  pairwise  harmonic  coupling  among  adjacent  oscillators)  failed  to  accurately 
reproduce  the  HCAO-simulated  spectral  regions. 

The  present  work,  which  was  carried  out  without  recourse  to  the  experimental  sarin  spectral 
regions,  illustrates  that  it  is  currently  feasible  to  predict  the  absorption  spectra  of  species  which 
are  difficult  to  synthesize,  handle,  or  otherwise  acquire.  In  addition  to  their  utility  in  guiding  ex¬ 
perimental  investigations,  the  simulated  overtone  spectral  regions  will  be  necessary  to  correctly 
assign  experimental  overtone  spectra  owing  to  the  large  number  of  similar  but  nonequivalent  CH 
oscillators  present  in  sarin, 
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local  mode,  overtone  transition,  Fermi  resonance,  coupling,  seating  factor,  normal  mode,  ab 
initio,  calculated  frequencies,  vibrational  mode,  molecular  vibration 
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